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Abstract. 7-ray binaries are colliding wind binaries (CWB) composed of a massive star a non- 
accreting pulsar with a highly relativistic wind. Particle acceleration at the shocks results in emis- 
sion going from extended radio emission to the 7-ray band. The interaction region is expected to 
show common features with stellar CWB. Performing numerical simulations with the hydrodynam- 
ical code RAMSES, we focus on their structure and stability and find that the Kelvin-Helmholtz 
instability (KHI) can lead to important mixing between the winds and destroy the large scale spiral 
structure. To investigate the impact of the relativistic nature of the pulsar wind, we extend RAMSES 
to relativistic hydrodynamics (RHD). Preliminary simulations of the interaction between a pulsar 
wind and a stellar wind show important similarities with stellar colliding winds with small relativis- 
tic corrections. 
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1. INTRODUCTION 

7-ray binaries share a common structure with coUiding wind binaries composed of two 
massive stars. In the latter, the interaction of the winds creates a double shock structure 
which geometry depends on the momentum flux ratio of the winds. Extensive numerical 
studies have studied the different instabilities that arise in the colliding wind region 
[1, 2]. At larger scale, a spiral structure is expected but its exact geometry is still 
being studied [3]. Relativistic simulations of 7-ray binaries reveal a structure similar 
to colliding stellar winds [4] and the development of instabilities [5]. Still, the exact 
impact of the relativistic nature of the pulsar wind has never been established. The aim 
of our work is to highlight the similarities and differences between 7-ray binaries and 
colliding stellar winds using numerical simulations. 

2. STELLAR COLLIDING WIND BINARIES 

RAMSES is a second order Godunov method, that solves the equations of hydrody- 
namics. We use adaptive Mesh Refinement (AMR) to increase resolution and properly 
model the instabilities while simulating the large scale structure at reasonable compu- 
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tational cost. The winds are generated following the method by Lemaster et al. [6]. We 
include a passive scalar to determine mixing between the winds. 
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FIGURE 1. Density and mixing for 2D simulations. In the first panels vi = 20v2, in the last panels 
vi = 1/20 = 0.05v2. The length scale is the binary separation. 

Fig. 1 shows density and mixing for winds with moment flux ratio M2V2/^iVi = 
0.0625 and inverted velocity ratios. The spiral structure has different properties accord- 
ing to the velocity of the dominant wind v\. Using additional simulations, we determine 
that the step of the spiral is mostly set by vi x P where P is the orbital period of the 
system. Still, the weaker wind can account for a non-negligible deviation to this value. 
There is a distinction between both spiral arms, clearly visible when vi = 20v2. The 
spiral arm propagating into the faster, lower density wind expands, while the arm prop- 
agating into the denser and slower wind gets compressed. Mixing is more important in 
the wider arm [7]. When vi/v2 = 20, the spiral structure is destroyed due to the KHI 
while the structure is stable when vi/v2 = 0.05 owing to the important density gradient 
between the winds. 

These simulations indicate the large scale structure of CWB is strongly dependant on 
the wind properties and that the KHI can account for the destruction of the large scale 
structure in some cases. It also leads to important mixing between the winds. In y-ray 
binaries, it could enhance cooling of non-thermal particles and affect the large scale 
synchrotron emission [8]. To determine the impact of the relativistic pulsar wind on the 
structure and stability of the colliding wind region, we have extended RAMSES to RHD. 



3. RELATIVISTIC HYDRODYNAMICS WITH RAMSES 

The equations of RHD can be written as a system of conservation equations (c = 1) 

dt^ dxj - ^ ( D \ ( Tp 

SM,^dm^ = With M U T^phvi 1 (1) 



dt ^ dx 
dt dxj - ^ 

where D is the mass density, M the momentum density and E the energy density in 
the frame of the laboratory. The subscripts ij stand for the dimensions, 5ij is the 
Kronecker symbol, h is the specific enthalpy, p is the proper mass density, v/ is the 
fluid three- velocity, P is the gas pressure and 7 the adiabatic index. The Lorentz factor F 
is given by F = (Vl —v^)~^/^. These equations have a similar structure to the equations 
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FIGURE 2. Left panel : Sod test = 1.8). Right panel : Simulation of the propagation of a 3D 
relativistic jet (Tmax = 7.1). From top to bottom: density at ^ = 20,30,40 in a 3D jet starting form the 
left boundary of the domain. 
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FIGURE 3. Left panel : density map of a simulation with equal momentum flux, with = 0.5, the star 
is on the left, the pulsar on the right. Right panel : position of both shocks and the contact discontinuity, in 
simulations with different values for the velocity of the pulsar wind. We have Vp= 0.01 (thin dotted line), 
0.1 (thick dotted lined), 0.5 (thick dashed line) and 0.9 (soHd line). 



of hydrodynamics but are more complex to solve because strongly coupled to each other 
by the Lorentz factor and the enthalpy. An additional numerical constraint arises from 
the the fact that the velocity must remain subluminal. The similarity with the equations 
of hydrodynamics allows us to closely follow the algorithm implemented in RAMSES, 
performing localised changes. 

We adapted the transition from the conservative variables (D,M,£')^ to the primitive 
variables (p, v,P)^ [9]. Second-order precision is implemented in RAMSES following a 
MUSCL-Hancock method and requires the determination of the Jacobian matrix of the 
system given by Eq. 1. The relativistic summation of velocities changes the determina- 
tion of the wavespeeds and the timestep. The implementation of RHD within the AMR 
structure requires adaptations when determining variables at a given refinement level 
/ by using the variables at level / — 1 or / + 1 . Our current implementation passes the 
standard numerical test. Fig. 2 shows the results of a 3D simulation of an axisymmetric 
jet following the setup by Del Zanna and Bucciantini [10]. This tests show satisfactory 
results and indicates the code is ready for scientific use. 



4. SIMULATING 7-RAY BINARIES 

Our aim is to model 7-ray binaries. We focus on the small scale structure of the interac- 
tion between a stellar wind and a pulsar wind. The goal is to understand the impact of 
relativistic effects both on the structure and stability of the interaction region.We neglect 
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orbital motion and focus on winds with equal moment fluxes. We perform preliminary 
simulations with various values of the momentum flux ratio and pulsar wind velocity. 
They prepare a large scale simulation to determine whether a stable structure is possible 
[5]. We also want to determine the Lorentz factor downstream, as it may account for 
boosted emission. Fig. 3 shows the density map for a simulation with the speed of the 
pulsar wind Vp = 0.5. The KHI develops in a similar fashion than in the classical case. 
The right panel shows the positions of the discontinuities for simulations with increas- 
ing values for the pulsar wind speed. The higher the value, the more the shocks are bent 
towards the star. This is a relativistic effect due to the impact of transverse velocities on 
the structure of shock 

5. PERPSECTIVES 

We performed high resolution simulations of colliding wind binaries at a spatial scale 
never reached before. We showed that the KHI may destroy the expected large scale 
structure. Simulations of WR 104 match well with the observed structure and indicate 
cooling has to be taken into account to allow dust formation in this system. To model 
7-ray binaries, we extended RAMSES to relativistic hydrodynamics. Preliminary sim- 
ulations of 7-ray binaries confirm a similar structure to stellar binaries. The relativistic 
extension of RAMSES allows the use of AMR and is suited for the study of gamma-ray 
bursts, relativistic jets or pulsar wind nebulae. It will be part of the next public release. 
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